Polynomial invariants and Bell inequalities 
as entanglement measure of 4-qubit states 
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We compare the polynomial invariants for four qubits introduced by Luque and Thibon, PRA 
67, 042303 (2003), with optimized Bell inequalities and a combination of two qubit concurrences. 
It is shown for various parameter dependent states from different SLOCC classes that it is possible 
to measure a genuine 4-qubit entanglement with these polynomials. 



PACS numbers: 03.67.Mn,03.65.Ud 
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The classification of multiqubit entanglement is actu- 
ally a wideley discussed field. Until now no measure for 
genuine n-qubit entanglement (n > 4) is known, in con- 
trast to the 2-qubit concurrence or the 3-qubit tangle. In 
this paper we will show by example that the combination 
of optimized Bell inequalities and polynomial invariants, 
introduced by Luque and Thibon 4], yields an entangle- 
ment measure for four qubits. 

In we argued that the comparison between an opti- 
mization of Bell-type inequalities [3| and a combination 
of the global entanglement measure Q [3 with the sum 
of the squared 2-qubit concurrences Gy ul yields the 
same parameter dependence of the entanglement mea- 
sures. We will show that one of the invariants shows the 
same behavior. 

The states which are in the focus of our research belong 
to different entanglement classes. Miyake and Ver- 
straete et al. discussed the SLOCC classification of 
multiqubit states. Miyake connected the classification to 
hyperdeterminants and showed that the representative of 
the outermost 4-qubit entanglement class is the state 

\G afi ys) = a(|0000) + |1111)) + /3(|0011) + |1100)) 

+ 7QOIOI) + |1010» +£(|0110) + |1001» (1) 

This state is equivalent to the state G a bcd of Vertraetes 
nine SLOCC families for four qubits. A criterion for this 
class is a nonzero hyperdeterminant DetA^, called A in 
the paper by Luque and Thibon. 

For the calculation of the SLOCC class affiliation, we 
take the easy-to-use criteria recently introduced by Li et 
al. 0- There starting points were representative states 
like the GHZ- or the W-state and an existence criterion 
for the SLOCC transformation. Two states \4>), are 
invariant under SLOCC if local invertible 2x2 matrices 
A, B, C, D exist with \<j>) = A® B ® C ® D\ip). Taking 
the state under consideration as \<j>) and the final state 
as one can easily calculate these criteria, that means 
equation under the condition that the matrices exist. 
In the following we discuss the results for three different 



classes of states. 

1) In our first example we will consider the 4-qubit pa- 
rameter dependent GHZ state, 



7|0000) + a/1-7 2 |1111) 



(2) 



with 7 e]0, 1[. It is assumed that this state shows 
genuine 4-qubit entanglement. The concurrences be- 
tween any pair of qubits, CV,-, are 0. They are calcu- 
lated the usual way yj |(j . The concurrence is defined as 
the maximum dj = maxj-y/A]" — \f\^, — VX? — V-^b 0}, 
and the \ are the eigenvalues of the matrix pij — 
Pij(cr v ® a y )p*j(a y ® a v ), with pij, the reduced density 
matrix to the qubits i and j. 

The global entanglement is Q = 47 2 (1 — -f 2 ). The cal- 
culation of the Luque/Thibon invariants yields H = 

0, D xt = 0, S = 7 4 (- 



7VT^, L 
7 2 ) 2 /12, T = 



M = N 



-1 



7 b (-l + 7 2 )7216 and A = (see ap- 



pendix). It is nicely seen that the 5* resp. T invariant 
and the global entanglement measure are the same, up 
to a square root and a constant factor. In Fig^we show 
the Global Entanglement Q, the invariant S and the Bell 
optimization as function of 7 2 . 

We use the Mermin-Klyshko-type Bell inequalities. You 
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FIG. 1: Optimized Bell inequality (black dots) on the right 
y-axis, Global Entanglement Q (dashed line) and S invariant 
(black line) on the left y-axis, as function of 7 2 for the 4-qubit 
parameter dependent GHZ state. The dotted lines show the 
Bell inequality condition for 3-qubit entanglement (> 8) resp. 
4-qubit entanglement (> 16). 
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get the corresponding operators out of the following re- 
cursion relations: 

F N = l -{D N + D' N )F N - 1 + l -{D N -D' N )F l N _ l (3) 
F' N = \(D N + D'M^ + \(D' N - D N )F N ^ (4) 

with F 2 = (A'B + AB') + (AB - A'B') and 
F-2 = {A'B + AB') - (AB - A'B'). The a''\B® 
resp. can be written as sums of Pauli matrices, e.g. 
= aW • <j* A for qubit A, with normalized vectors a^'\ 

These Fji operators are optimized and yield a criterion 
for 4-qubit entanglement, if (F 4 ) 2 + (i^) 2 > 16 HQ. 
Though the Bell optimization for the parameter depen- 
dent GHZ state has some difficulties near 7^0 and 
7 ~ 1 as described in |2(, the three quantities match 
nicely, especially quantified at the maximum 7 2 = 1/2. 
The state does not belong to the outermost 4-qubit 
SLOCC class as described by Miyake, because the 
hyperdeterminant criterion is not fulfilled, A = 0. 
The simple criteria for the GHZ class after Li et al.is 
-7-v/l - 7 2 and is fullfilled for 7 e]0, 1[. 

2) The classification of the next state is more com- 
plex. Parts of it were already done in Q. The state is 
one of the eigenstates of a special 4-qubit Heisenberg 
model, with antiferromagnetic coupling constants J and 
J s bot > 0, and is written as: 

\<h) = ft (-I0011) + |oiio) - |iooi) + |iioo)) 

+ /3 2 (-|0101) + 11010)) (5) 

with the two amplitudes fa and fa given by fa = 
(4 + (-J + 2J S + 5) 2 /(2J 2 )Y 1 '\ fa = (4 + (-J + 

2J S + 5) 2 /(2J 2 )) 1/2 J/{2S) and the abbreviation 6 = 
(9J 2 -4JJ s + 4J 2 ) 1 /2. 

We know from the optimization of Bell-type inequali- 
ties that this state fullfills the criteria for 3- resp. 4-qubit 
entanglement. The global entanglement measure is con- 
stant, Q = 1. Additionally the parameter dependence of 
the optimized inequality is similar to the course of the 
combination Q — C 2 j = 1 — ^ C-L The parameter 
dependence of both is shown in Fig. [z] as a function of 
J with constant J s — 2 resp. of J s with constant J = 2. 
The calculation of the Luque/Thibon invariants yields 
the following. H = -2/3 2 - (3 2 , L = -N = 0{(0{ - /3|), 
M = and D xt — ~fa\fi\ (see appendix). The more 
interesting ones are 

S=^(-^l + P 2 2) 2 (6) 

T= 2l^ 6 (-4/?? +/3 2 2 ) 3 (7) 

and A = 0, with the general relation A = S 3 - 27T 2 . 
In Fig. |21 we also show the paramater dependence of the 
invariant S, which is normalized to 1. All three quanti- 
ties plotted show the same behavior. The coincidence 




FIG. 2: Optimized Bell inequality (line with dots) on the 
right y-axis, 1 ~Y1 Cij (dashed line) and S invariant (black 
line) as function of J resp. J s for the state [^2) on the left 
y-axis. The dotted line shows the Bell inequality condition 
for 4-qubit entanglement. J and Js are complicated functions 
of Pi and 02 (see text). 



is very nice seen at the maxima resp. minima. For 
J = 0, J s = 2 resp. J — 2,J S ^ 00 all three quan- 
tities show a maximal 4-qubit entanglement. Also the 
minimum at J = 2, J s = 2 matches very nicely. 
Since the invariant A is equal to 0, the state docs 
not belong to the Gup^s class. But the state fullfills 
for almost all paramters, except a small range around 
J = J s = 2, the Bell condition for 4-qubit entanglement. 
Because of this fact we will test the classification criteria 
for the GHZ SLOCC class as introduced by Li et al. [jj. 
We get the following equations: 

2/3 2 +/3 2 V0, -Pf=0 and ft ft = 0. (8) 

These are solved for fa = and fa 7^ 0. This 
is the case for the limits J s = 2, J — > and 
J = 2, J s — > 00. The state is reduced to a GHZ 
type state \fa) = (— 10101) + |1010))/V2, and belongs to 
the GHZ SLOCC class. In all other cases, that means 
fa 7^ 0, the state belongs to another not yet quantified 
SLOCC class. 

3) In the next part we study a state with A 7^ 0. 
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FIG. 3: Optimized Bell inequality (line with dots) on the 
right y-axis, 1 ~Y1 @ii (dashed line) and S invariant (black 
line) as function of 7 for the state \G a -y), on the left y-axis. 
The dotted line shows the Bell inequality condition for 4-qubit 
entanglement. 



Therefore we take the Miyake state \G a pys) an d reduce 
the number of paramters. We set (3 = S = 7 and 
2a 2 + 67 2 = 1: 



\G ai ) = a(|0000) + |1111» +7(|0011) + 11100)4 
|0101) + |1010) + |0110) + |1001>), 



(9) 



and we choose a resp. 7 real, 7 e [0,1/v^]. A gain we cal- 
culate the Luque/Thibon invariants and compare them 
with the Bell optimization and the combination of global 
entanglement and squared concurrences. The global en- 
tanglement is constant Q = 1 and the concurrences are 
all equal, C\% — C13 = C14 = C 23 = C24 = C34, with 



C12 — 



A"f{a — 7) 
2( 7 2 -a 2 ) 



if 
if 



2V2 ^ 7 - V6 



(10) 



with a defined above. 

For the Luque/Thibon invariants we get the following: 
H = 1/2, L = M = N = 0, A,* = 1/4 7 2 (1 - 8 7 2 ) 2 and 



S= — -i 7 2 (i- 8 7 2 ) 2 
192 4 ' v 1 ' 

1 7 2 7 7 4 

T = — + — - 

13824 192 48 



24 7 8 - I287 
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(6 7 2 -l)(24 7 2 -l) 2 (87 3 - 7) 6 



(11) 

(12) 
(13) 
(14) 



In Fig. |3|we show the parameter dependence of our mea- 
sures. It is again clearly seen that the course of the in- 
variant S, which is again normalized to 1, the optimized 
inequality and the function of the squared concurrences 
yield the same result in the parameter dependence. Es- 
pecially for 7 = 0, the state reduces to a GHZ state, and 



all three quantities have a maximum. Also the minimum 
at 7= 1/2 and the maximum at 7 = l/2y6 match. 
The state \G a pjs) is the representative of the outermost 
SLOCC entanglement class with the criteria A 7^ 0. 
Also for our special choosen state \G ai ) the hyperde- 
tcrminant A is different from 0, except for three values, 
7 = ±l/V6,±l/\/8,±l/\/24. We now take the criteria 
for the GHZ SLOCC class by Li et al. 



3 7 2 ^ A 



2 2 
a 7 



7 







(15) 



and test them with the calculated roots for A. For 
7 = l/v§ these equations are solved and the state 
belongs to the GHZ class. 

Another interesting feature of the invariants is found 
from the comparison of the invariant H and the inequal- 
ity which comes out of the criteria for the GHZ SLOCC 
class. Up to a minus sign, they are equal in all our 
examples. That means, a nonzero invariant H is a crite- 
rion for the affiliation of a state to the GHZ SLOCC class. 

Conclusions and discussions - It is shown in this 
paper that the genuine 4-qubit entanglement could be 
measured with a polynomial invariant. This invariant 
S yields the same parameter dependence as optimized 
Bell-type inequalities and a combination of global 
entanglement and 2-qubit concurrences for the states we 
have choosen from different SLOCC classes. 



APPENDIX: LUQUE/THIBON INVARIANTS 

Here we give the general equation for the 
Luque/Thibon invariant D xt , beacause it is not ex- 
plicitly calculated in [y. A pure 4-qubit state can be 
written in the computational basis as = X)iio a *N)' 
The invariants can then be expressed in terms of the ai. 
Here especially the invariant D xt : 



Act = (-a 11 a 13 +a 15 a 9 )(-(a 3 a 4 + a 2 a 5 - a x a 6 - a a 7 ) 
(— a ai4 + ai 2 a 2 + ai a4 — a 6 a s ) + (a 2 a 4 — a a 6 ) 
(-aia 14 - a ai 5 + ai 3 a 2 + ai 2 a 3 + 

ana4 + a 10 a 5 - a 7 a 8 - a 6 a 9 ))+ 
(~a w ai 2 + ai 4 a s )(-(a 3 a 5 - aia 7 ) 
(— aiai4 - a ai 5 + a 13 a 2 + ai 2 a 3 + 
011O4 + aloas - a 7 a 8 - a 6 a 9 )+ 
(a 3 a4 + a 2 a 5 — aia 6 — aoa 7 ) 
(-aids + ai 3 a 3 4- ana 5 - a 7 a 9 ))- 
(-011012 - a w ai 3 + ai 5 a 8 + a u a 9 )((a 3 a 5 - 0,10,7) 
(a ai4 — a\ 2 a 2 — ai a 4 + a 6 a 8 ) + (— a 2 a 4 + a a 6 ) 

(01O15 - ai 3 a 3 — on 05 + 07(19)) (A.l) 
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